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ABSTRACT 

Aims. To study the high-mass star-forming process, we have started a large project to unveil the gas kinematics close to young 
stellar objects (YSOs) through the Very Long Baseline Interferometry (VLBI) of maser associations. By comparing the high spatial 
resolution maser data, tracing the inner kinematics of the (proto)stellar cocoon, with interferometric thermal data, tracing the large- 
scale environment of the hot molecular core (HMC) harbouring the (proto)stars, we can investigate the nature and identify the sources 
of large-scale motions. The present paper focuses on the high-mass star-forming region G16. 59-0.05. 

Methods. Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser 
species in G16.59-0.05: H z O at 22.2 GHz (4 epochs), CH 3 OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, 
we performed high-resolution (> 071), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the 
star-forming region at 1.3 and 3.6 cm. 

Results. This is the first work to report accurate measurements of the relative proper motions of the 6.7 GHz CH3OH masers. The 
different spatial and 3-D velocity distribution clearly indicate that the 22 GHz water and 6.7 GHz methanol masers are tracing different 
kinematic environments. The bipolar distribution of 6.7 GHz maser l.o.s. velocities and the regular pattern of observed proper motions 
suggest that these masers are tracing rotation around a central mass of about 35 M G . The flattened spatial distribution of the 6.7 GHz 
masers, oriented NW-SE, suggests that they can originate in a disk/toroid rotating around the massive YSO which drives the 12 CO (2- 
1) outflow, oriented NE-SW, observed on arcsec scale. The extended, radio continuum source observed close to the 6.7 GHz masers 
could be excited by a wide-angle wind emitted from the YSO associated with the methanol masers, and such a wind is proven to 
be sufficiently energetic to drive the NE-SW 12 CO (2-1) outflow. The H 2 masers distribute across a region offset about 075 to the 
NW of the CH3OH masers, in the same area where emission of high-density molecular tracers, typical of HMCs, was detected. We 
postulate that a distinct YSO, possibly in an earlier evolutionary phase than that exciting the methanol masers, is responsible for the 
excitation of the water masers and the HMC molecular lines. 

Key words. Masers - Techniques: high angular resolution - ISM: kinematics and dynamics - Stars: formation - Stars: individual: 
IRAS 18182-1433, G16.59-0.05 
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~j \ 1 . Introduction Baseline Interferometry (VLBI) technique to determine the po- 

3t sition of the maser spots with milli-arcsec accuracy, and com- 

O . The process of massive star formation is based on a fine tun- ^ positions at different VLBI epochs gives accurate mea- 

7H ■ ing between the gravitational force, that sets in the collapse surement of the spot proper motions . Combining proper mo- 

> , of a Jeans-critical cloud, against the different types of pres- tions whh ±e maser line . of . sight (Lo s) velocit ie S (derived via 

sures, thermal, magnetic, turbulence, and radiation, that regu- Doppler shift of the observed maser freqU ency), one can recon- 

late the time scale of the accretion (for recent reviews see, e.g., stmct the m 3 _ D kinematics of the masing gas . with this in 

b . | Zinnecker&Yorke| | 2007t | Beuther et al. || 2007|) . To characterize mind> we haye started a large project whose final aim . g draw . 

the relationship between accretion-ejection phenomena it is es- ing a compre hensive picture of the dynamical processes in high- 

sential to resolve the dynamical structures close to the central mass star . forming regions (HMSFRs), by relating the pc-scale 

engine. Maser emission from several molecular lines is an use- phenomenai traced by thermal continuum and mo i eC ular lines, 

ful signpost of the hot, dusty environment (i.e. hot molecular with the AU . scale kine matics of the gas around young stellar 

core, HMC) where massive star formation takes place, which, objects (YSOsX traced by maser emissions . 
because of the extremely high dust extinction, could not be 

investigated with optical and near-infrared (NIR) diagnostics. In particular, observing targets where different molecular 

Maser emission is concentrated in narrow (usually < 1 km s _1 maser species appear to be spatially associated, we can hope to 

broad), strong (up to 10 6 Jy) lines and arises from compact emis- better sample the (proto)stell ar environment (e.g., iGoddi et aTl 

sion cente rs ("maser spots") which have typical si zes of a few I2007t iMoscadelli et al.ll2007t) . As a first step, we focus on the 

AU (e.g JMiniere"tail2 002: Moscadelli et al. 2003). The strong three most powerful maser transitions, that of water (H2O) at 

brightness of the maser emission allows us to use the Very Long 22.2 GHz, methanol (CH3OH) at 6.7 GHz, and hydroxyl (OH) 

at 1.665 GHz. H2O masers typically trace fast shocks (of sev- 

Send offprint requests to: A. Sanna, e-mail: asanna@ca . astro . it eral tens of km s _1 ) excited by either collimated jets found at 
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the base of molecula r outflows, or wide-angle winds emitted 
from the YSOs (e.g. iGwinn et al.lll99l iTorrelles etai1l2003t 
iMoscadelli et al.l2007l) . OH masers are likely associated with the 
slow expandin g (typically a few km s" 1 ) ionization front of H n 
regions (e.g.. lFish& Reidll2007h . The origin of CH3OH masers, 
instead, is currently a matter of debate. To date, single-epoch 
VLBI observations have provided accurate spatial and l.o.s. ve- 
locity distribution of the 6.7 GHz masers towards a few tens of 
sources. To account for the ordered gradient of l.o.s. velocities 
observed in some objects, three hypotheses have be en proposed: 
circumstellar (Keplerian ) discs seen ed ge-on (e.g., Norris et al.l 
U998tlMinier et al.l2000l) ; outflows (e.g. JDe Buizerl2003l) : prop- 
agating shock f ronts through rotating de nse cores and ring-like 
structures (e.g.. [Bartkiewicz et al.ll2009h . The only way to dis- 
criminate among these different interpretations is deriving maser 
proper motions via multi-epoch VLBI observations. 

The present paper focuses on the HMSFR G16.59-0.05. In 
Sect. 2, we provide an up-to-date review of the single-dish and 
interferometric observations towards this region. Section 3 de- 
scribes our VLBI observations of the 22.2 GHz H2O, 6.7 GHz 
CH3OH, and 1.665 GHz OH masers, together with the new Very 
Large Array (VLA) A-configuration observations of the radio 
continuum emission at 1.3 and 3.6 cm, which complemented 
VLA-C archival data at 0.7, 1.3, and 3.6 cm. Details of our data 
analysis are given in Sect. 4. In Sect. 5, we illustrate the spa- 
tial morphology, kinematics, and time-variability of individual 
maser species, and present results from our and archival VLA 
observations, constraining the properties of the radio continuum 
observed in correspondence of the masers. Section 6 discusses 
the spatial association of the maser species and their overall 
kinematics, and draws a comprehensive picture of the phenom- 
ena observed in the HMSFR G16.59-0.05 on angular scales 
from a few mas to tens of arcsec. The main conclusions are sum- 
marized in Sect. 7. 



2. The HMSFR G1 6.59-0.05 

The HMSFR G16.59-0.05 (IRAS 18182-1433) has a near/far 
kinematic distanc43 of 4.4/11.9 kpc evaluated assuming a sys- 
temic velocity (V^) with respect to the local standard of 
rest (LSR) of 59.9 km s _1 , inferred fr om observations of 
the NH 3 (3 3) and CH 3 CN (5-4) lines dCodella et all \l99% 
IFuruva et alj [2008). For consistency with previous analyses 
of this source, in the following the near kinematic distance 
is adopted. At the near kinematic distance, the bolometric 
luminosity calculated from the IRAS fluxes is 2 x 10 4 L Q 
( Sridhara n et al.l 120021) and th e large-scale clump mass is > 
10 3 Mm dBeuther et all I2002at iFaundez et all l2004t iHill et alj 



120051: IFuruva et al.ll2008l)T 

On a large scale, mid-IR observations toward G16.59-0.05 
resolved the IRAS source in the region detecting 2 objects sepa- 
rated by about 10" along the NW-SE direction dDe Bu izer et al. 
120051: IFuruva et al.ll2008l Fig. 1). While the SE emis sion is asso- 
ciated with a MSX point source (IWalsh et al] l2003). the fainter 
NW source 0.2 J y at 18.1^m) is a ssociated with weak, red- 
dened NIR emission dTesti et al.ll 19941) and wa s tentatively clas- 
sified as a Class I object ( De B uizer et al. 20051). The NW source 
is also coincident with the submillimeter (450 and 850 /urn) con- 
tinuum emission from the region, which may indicate that this 
is a younger, more embedde d source than the one responsible 
for the SE emission (Wals h et alj 120031: IWilliams et alj [2004. 



120051) . Single-dish surveys toward G16.59-0.05 detected red- 
shifted self-absorption in N 2 H + , HCO + , H 13 CO + , and H 2 CO ro- 
tational line s , which suggests the p r esence of inflow motion s 
dFulleretal.1 120051: IWu et alJ 120071: [Thomas & Fuller! 12008). 
Bolometer array ob servations at 1.2 mm (Beu ther et al.l l2002a: 
Faun dez et al.| [2004) show a peak of dust emission in correspon- 
dence of the two mid-IR sources and an additional fainter source 
about 100" toward the S E, identified as a candid ate High Mass 
Starless Core (HMSC) bv lSridharan et all d2005l) . 

High resolution observations at millimeter and radio wave- 
lengths, both in line and continuum emission, revealed multiple 
sources within the mid-IR region (see Fig. [8] where we sum- 
marize both previous observations and our results). A molecu- 
lar core with a kineti c temperature of 54 K is identified through 
NH 3 observations bv lCodella et alJd!997l) . Sensitive VLA-C ob- 
servations at 3.6, 1.3 cm, and 7 mm detected two weak cen- 
timeter sources and one 7 mm co ntinuum source, labeled "c", 
"b", and "a" bv lZapatae t al. (2006), respectively. The centime- 
ter source "c" appears to be associated with the strongest mid-IR 
source; the other two VLA sources (with the 7 mm source "a" 
observed about 1" to the NW of the nearby centimeter source 
"b") are offset by about 10" to the NW of the strongest mid-IR 
source, and correspond in position with the NW mid-IR emis- 
sion. Continu um measurements a t 3 mm dFuruva et alj 1*2008) 
and 1.3 mm dBeuther et al.l 12006) peak also at the same posi- 
tion as the VLA source "a". The spectral index suggests that 
the millimeter emission traces the main massive proto-stellar ob- 
ject in the field, whereas the centimeter feature "b" might indi- 
cate a thermal jet o r, alternatively, an optically thin H 11 region 
(Zap ata et al.l 12006]) . Emission in several molecular lines typi- 
cal of HMCs was detected at an intermediate position bet ween 
source "a" and "b" (Beut her et al.l2 006; Furu va et al.l2 008). The 
rotational temperature inferred from the CH3CN lines is about 
130-150 K in this region, consistent with the presence of an 
HMC. 

In this work we will focus our attention on the NW sources, 
namely the HMC and the associated sources "a" and "b", hosting 
the maser emissions detected toward G16.59-0.05. OH masers 
were reported in the main 1 .665 and 1 .667 GHz, and the satellit e 
1.720 GHz lines (ICasweU & Havnesl 119831: lEdris et al.l [2007b . 
whereas the 1.612 GHz OH satellite was observed in absorption. 
CH3OH ther mal fines were observ ed at 25, 96.7, 157, 241.7, and 
255 GHz by iLeurini et alJ d2007l) , and maser emission was re- 
ported for the Class II 6 7 GH z dMentenl [19911: [ Caswell et al 
I l995bt ISzvmczak etail l2000h and 12 2 GHz dcaswell et al. 
1995aD, an d the Class I 44 1 GHz dSlvsh et al.l 11994 and 
95.2 GHz dVal'tts et al.l l2000t) transitions. High-resolution (> 
1"), interferom etric observations towards G 16.59-0.05 of th e 
22.2 GHz H 2 dForster & Caswel|l 19991: feeuther et al.l[2002d) . 
1.665 GHz OH dForster & Caswell 19991) . and 6.7 GHz CH 3 OH 
dWalshetal.1 19981) " maser lines, associate the maser activity with 
the sources "a" and "b". 



1 We adopt a "r evised" kinematic distance using 
iReid etal.U2009bl) 



A wide angle, massive (> 30 M Q ) bipolar outflow emerg- 
ing from the HMC was detected through 12 CO and CO isotopo- 
logues rotational lines (1-0 a nd 2-1) at an angular resolution of 
about 10" by several authors (|Beuther et all2002bt IFuruva et al.l 
2008; Lopez-Sepulcre et al. 2009). At the higher spatial resolu- 
tion (« 3") of the SMA interferometer, the 12 CO (2-1) emissio n 
was resolved in a multiple outflow system (Beut her et al.l l2006): 
the strongest and better collimated lobes are oriented along two, 
approximately perpendicular directions, NW-SE and NE-SW. 
the prescription of The outflow system is centered on the position of the millime- 
ter source. Furthermore, VLA SiO (1-0) observations revealed 
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emission elongated toward the N-S direct ion, which could be 
suggestive of an additional third outflow (Be uther et al.ll2.Q06h . 

The large collection of multi-wavelength data reported here 
suggests that the G16.59-0.05 region is an active site of multiple, 
massive star formation. 

3. Observations and Calibration 

3.1. Archival Data and VLA continuum 

We observed G16.59-0.05 with the VLA0in the C- and A-array 
configurations at X band, and in the A-array configuration at K 
band. The VLA-C and VLA-A observations were made respec- 
tively in April and October 2008. The source had already been 
observed at K and Q bands with the VLA-C by Zapata et al. 
(2006). We retrieved their data from the VLA archive (Program 
Code: AZ149) and re-reduced them to improve on the angular 
resolution of the maps presented in their paper. In particular, we 
imaged the Q-band data using natural weighting, achieving an 
angular resolution about a factor 2 better than that of the tapered 
maps by Zapata et al. (2006). At 3.6 cm we used the continuum 
mode of the correlator, resulting in an effective bandwidth of 
172 MHz. For the 1.3 cm observations, we used mode "4" of the 
correlator, with a pair of 3.125 MHz bandwidths (64 channels) 
centered on the strongest H2O maser line and a pair of 25 MHz 
bandwidths (8 channels) sufficiently offset from the maser lines 
to obtain a measurement of the continuum emission. This mode 
was used to improve on the relative position accuracy between 
the maser emission and the high-resolution observations of the 
1.3 cm continuum. The two bandwidths centered at the same 
frequency (measuring the two circular polarizations) were aver- 
aged. 

At X-band, 3C286 (5.2 Jy) and 3C48 (3.1 Jy) were used as 
primary flux calibrators, while 1832-105 (1.4 Jy) was the phase 
calibrator. For the K-band observations, the primary flux calibra- 
tor was 3C286 (2.5 Jy), the phase calibrator 1832-105 (1.0 Jy), 
and the bandpass calibrator 1733-130 (3.7 Jy). 

The data were calibrated with the NRAO AIPS software 
package using standard procedures. Only for the VLA-A data 
at 1.3 cm, several cycles of self-calibration were applied to the 
strongest maser channel, and the resulting phase and amplitude 
corrections were eventually transferred to all the other line chan- 
nels and to the K-band continuum data. This procedure resulted 
in a significant (at least a factor 2) improvement of the signal-to- 
noise ratio (SNR). 

3.2. Maser VLBI 

We conducted VLBI observations of the H2O and CH3OH 
masers (at several epochs), and of the OH maser (single epoch) 
toward G16.59-0.05 in the K, C, and L bands, respectively. To 
determine the maser absolute positions, we performed phase- 
referencing observations by fast switching between the maser 
source and the calibrator J1825-1718. This calibrator has an 
angular offset from the maser source of 2.9 and belongs to 
the list of sources defining the International Celestial Reference 
Frame (ICRF). Its absolute position is known to better than 
+2 mas and its flux measured with the Very Long Baseline Array 
(VLB A) at S and X bands i s 48 and 112 mJy beam 1 , respec- 
tively dFomalont et alJl2003l) . Five fringe finders (J1642+3948; 

2 The VLA is operated by the National Radio Astronomy 
Observatory (NRAO). The NRAO is a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



J1751+0939;J1800+3848;J2101+0341;J2253+1608) were ob- 
served for bandpass, single-band delay, and instrumental phase- 
offset calibration. 

Data were reduced in AIPS employing the VLBI spectral 
line procedures. Visibility amplitudes were calibrated applying 
the information on the system temperature and antenna sensi- 
tivity provided by the monitoring procedure at each VLBI an- 
tenna site. A priori delay and phase calibration for feed rotation, 
correction for ionospheric delays (AIPS task TECOR), and for 
inaccuracies of the extrapolated Earth's orientation parameters 
used at the time of correlation (AIPS task CLCOR), were ap- 
plied to the data set (e.g., EVN data analysis guid^3). A pos- 
teriori phase calibration for short-term atmospheric fluctuations 
and phase-reference structure was derived by fringe-fitting and 
self-calibrating the strong, compact, emission in a given maser 
channel. Applying this calibration, visibilities of all maser chan- 
nels were referred in phase (i.e. in position) to the emission cen- 
troid of the reference maser channel. The determination of the 
absolute position of the reference maser channel was accom- 
plished in either of the two following ways: imaging the calibra- 
tor J 1825- 17 18 data after applying corrections derived working 
with the reference maser channel (inverse phase-referencing); 
or, by performing the reverse procedure, imaging the reference 
maser channel by applying the fringe-fitting and self-calibration 
solutions from the calibrator data (direct phase-referencing). In 
the inverse phase-referencing procedure, before fringe-fitting the 
maser data, the visibility phase center was shifted (AIPS tasks 
CLCOR and UVFIX) close to the effective maser posi tion, in or- 
der to avoid degradation of the calibrator image (e.g., Rei d et all 
2009a). For well detected (maser and calibrator) signals (SNR > 
10), the two procedures always gave consistent results. 



3.2.1. VLBA observations: 22.2 GHz H 2 masers 

We observed the HMSFR G16.59-0.05 (tracking center: 
R.A.(J2000) = 18 ;i 21 m 09!ll and Dec.(J2000) = -14°31'48"5) 
with the VLBAQ in the 616 — 5 2 3 H 2 transition (rest frequency 
22.235079 GHz). The observations (program code: BM244) 
consisted of 4 epochs: on April 9, June 28, and September 
18, 2006, and on January 4, 2007. During a run of about 6 h 
per epoch, we recorded the dual circular polarization through a 
16 MHz bandwidth centered on a LSR velocity of 60.0 km s _1 . 
The data were processed with the VLBA FX correlator in 
Socorro (New Mexico) using an averaging time of 1 s and 1024 
spectral channels. The total-power spectrum of the 22.2 GHz 
masers toward G16.59-0.05 is shown in Fig.Q](top panel). This 
profile was obtained by averaging the total-power spectra of all 
VLBA antennas, weighting each spectrum with the antenna sys- 
tem temperature (T jyj ). 

The natural CLEAN beam was an elliptical Gaussian with a 
FWHM size of about 1.4 mas x 0.4 mas at a PA. of -16° (east 
of north), with little variations from epoch to epoch. The inter- 
ferometer instantaneous field of view was limited to about 2'.'7 . 
At the various observing epochs, using an on-source integration 
time of about 2.5 h, the effective rms noise level of the channel 
maps (cr) varied in the range 0.01-0.04 Jy beam -1 . The spectral 
resolution was 0.2 km s _1 . 



3 www.evlbi.org/user_giiide/guide/user_guide.html 

4 The VLBA is operated by the NRAO. 
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3.2.2. VLBA observations: 1.665 GHz OH masers 

We observed the HMSFR G16.59-0.05 (tracking center: 
R.A.(J2000) = 18 A 21 m 09 s .21 and Dec.(J2000) -14°31'48'.'3) 
with the VLBA in the 2 n 3/2 J = 3/2 OH transition (rest 
frequency 1.665401 GHz), on April 13, 2007 (program code: 
BM244M). During a run of about 6 h, we recorded the dual cir- 
cular polarization through two bandwidths of 1 MHz and 4 MHz, 
both centered on a LSR velocity of 60.0 km s _1 . The 4 MHz 
bandwidth was used to increase the SNR of the weak L-band 
signal of the continuum calibrator. The data were processed with 
the VLBA FX correlator in two correlation passes using either 
1024 or 512 spectral channels for the 1 MHz and 4 MHz bands, 
respectively. In each correlator pass, the data averaging time was 
2 s. The T S y S -weighted mean of antenna total-power spectra for 
the right and left circular polarizations at 1 .665 GHz are shown 
in Fig. Q] (bottom panel). 

The natural CLEAN beam was an elliptical Gaussian with a 
FWHM size of 24 mas x 10 mas at a RA. of 3°. The interferom- 
eter instantaneous field of view was limited to about 18'.'5. With 
an on-source integration time of about 1.9 h, the effective rms 
noise level on the channel maps was about 0.01 Jy beam 1 . The 
1 MHz band spectral resolution was 0.2 km s _1 . 



3.2.3. EVN observations: 6.7 GHz CH 3 OH masers 

We observed the HMSFR G 16.59-0.05 (tracking center: 
R.A.(J2000) = lS h 2rW.13 and Dec.(J2000) = -14°31'48'.'5) 
with the European VLBI Network (EVNfl in the 5i - 6 A + 
CH3OH transition (rest frequency 6.668519 GHz). This work 
is based on 3 epochs (program codes: EM061, EM069), sep- 
arated by about 1 yr, observed on February 26, 2006, on 
March 16, 2007, and on March 15, 2008. At the first two 
epochs, antennas involved in the observations were Cambridge, 
Jodrell2, Effelsberg, Hartebeesthoek, Medicina, Noto, Torun 
and Westerbork. Since the longest baselines involving the 
Hartebeesthoek antenna (e.g., the Ef-Hh baseline is about 
8042 km) heavily resolve the maser emission and do not produce 
fringe-fit solutions, the Hartebeesthoek antenna was replaced 
with the Onsala antenna in the third epoch. During a run of about 
6 h per epoch, we recorded the dual circular polarization through 
two bandwidths of 2 MHz and 16 MHz, both centered on a LSR 
velocity of 60.0 km s _1 . The 16 MHz bandwidth was used to 
increase the SNR of the weak continuum calibrator. The data 
were processed with the MKIV correlator at the Joint Institute 
for VLBI in Europe (JIVE - Dwingeloo, The Netherlands) using 
an averaging time of 1 s and 1024 spectral channels for each 
observing bandwidth. The Effelsberg total-power spectrum at 
6.7 GHz toward G16.59-0.05 is shown in Fig. Q] (middle panel). 

The natural CLEAN beam was an elliptical Gaussian with 
a FWHM size of about 15 mas X 5 mas at a RA. of 30°, with 
little variations from epoch to epoch. The interferometer instan- 
taneous field of view was limited to about 9'.'2. At each observing 
epoch, using an on-source integration time of about 2.2 h, the ef- 
fective rms noise level of the channel maps varied in the range 
0.008-0.1 Jy beam 1 . The 2 MHz band spectral resolution was 
0.09 km s" 1 . 



5 The European VLBI Network is a joint facility of European, 
Chinese and other radio astronomy institutes founded by their national 
research councils. 



4. Data Analysis 

Mapped maser channels were searched for emission above a 
conservative threshold, equal to the absolute value of the min- 
imum in the map, typically corresponding to values greater than 
5-7 cr. Parameters (position, intensity, flux, and size) of the de- 
tected maser (and calibrator) emission have been derived by fit- 
ting a two-dimensional elliptical Gaussian. The term "spot" is 
used to refer to maser emission on a single channel map. We 
checked the contour plot of each detected spot and repeated the 
Gaussian fit to optimize the number of Gaussian components 
in case of complex spatial structure. We took particular care to 
identify and remove spurious spots owing to side-lobe effects of 
the synthesized beam. The uncertainty of the relative fit position 
is estimated using the expression: (A6) 2 = (0.5 FWHM/SNR) 2 + 
(50//as) 2 . The first contri bution represents the Gaussian fit un- 
certainty (iReid et al.lll988l) . with FWHM being the (undecon- 
volved) full-width half-maximum size of the spot, and SNR the 
signal-to-noise ratio of the Gaussian fit. The second contribu- 
tion represents a noise level added to take into account position 
uncertainties associated with cor relator-unmodeled signal prop- 
agation through the troposphere (Pradel et al. 2006|). 

The VLBI angular and velocity resolutions are sufficiently 
high to resolve the (spatial and velocity) structure of sin- 
gle masing-clouds. We collected spots observed on contiguous 
channel-maps and spatially overlapping (within their FWHM 
size), into a single maser "feature". The position and LSR ve- 
locity of a given feature are estimated from the error-weighted 
mean position and the intensity-weighted mean LSR velocity of 
its spots, respectively. If, for a given feature, spot position er- 
rors are comparable with the size of the spot distribution (as it 
is the case of most 22.2 GHz H2O masers, extended for a few 
tenth of mas), the feature's position uncertainty is evaluated by 
the error-weighted standard deviation of the spots positions. For 
extended maser features (as most 6.7 GHz CH3OH masers, ex- 
tended for a few mas), for which the diameter of the spot distri- 
bution is significantly larger than the spot position uncertainties, 
the position error is taken equal to the error-weighted mean of 
the spot position errors. The uncertainty in the absolute positions 
of the maser features is calculated by taking the square mean of 
three independent error terms: 1) the feature's relative-position 
uncertainty; 2) the position error of the (calibrator or reference 
maser) signal in the phase-referenced maps; 3) the position un- 
certainty of the calibrator J1825— 1718. This latter term domi- 
nates the absolute-position error for water and methanol maser 
features. 

We have used two main criteria to establish the correspon- 
dence of maser features over the epochs: 1) persistence of the 
relative distribution of a group of features (see inset in Fig. [3j 
see also God di et al.ll2006l); 2) assumptio n of uniform motions 
(e.g. jReid et al.lll988l;lGwinn et al.lll992l) . The tolerance in the 
change of peak LSR velocity for persistent features was fixed to 
less than their spectral FWHM (see Fig. [6}. Relative proper mo- 
tions for features persisting over (at least) 3 epochs, have been 
calculated by performing a linear fit of their varying (relative) 
position with time. The derived proper motions are a measure of 
the feature mean motion over a time baseline of (up to) 1 yr, for 
water, and 2 yr, for methanol masers. 

Absolute velocities are derived by adding to the relative ve- 
locities the absolute motion of the reference maser feature (the 
one to which velocities have been referred), determined with re- 
spect to the calibrator J 1825- 17 18. To estimate maser motions 
relative to the LSR reference frame of the HMSFR under study, 
the measured absolute motion of the reference maser feature has 
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to be corrected for the apparent proper motion due to the earth 
revolution around the Sun (parallax), the Solar Motion and the 
differential Galactic Rotation between our LSR and that of the 
HMSFR. Taking the IAU values for the Solar Motion and the 
Galactic Rotation (Ro = 8.5 kpc, ©o = 220 km s _I ), and assum- 
ing the near kinematic distance of 4.4 kpc, we derive an appar- 
ent motion of 12 km s to the west and 38 km s _I to the south. 
These values can be affected by large errors of 10-20 km s _I , 
owing to uncertainties in the Solar and Galactic standard values, 
as well as to the recent observational evidence that HMSFRs 
might be orbiting the G alaxy about 15 km s _1 slower than ex- 
pected for circular orbits (Rei d et alj200 9b). Absolute velocities 
have been estimated for water masers only, whose measured rel- 
ative motions (see Sect. 15.21 ) are found on average significantly 
larger than the expected uncertainty of the correction for the ap- 
parent motion. 



5. Results 

5.1. Radio continuum emission 

In the followin g, we focus our atten tion on the centimeter source 
labeled "b" bv lZaoata et all J2006l) . the one located nearby the 
millimeter continuum, since it appears to be spatially associ- 
ated with the observed 6.7 GHz CH 3 OH and 22.2 GHz H 2 
maser distributions in G16. 59-0.05 (Fig. [3]). Our VLA-A and 
VLA-C measurements of the 1.3 and 3.6 cm continuum, and 
the VLA-C archival data at 0.7 and 1.3 cm, are summarized in 
TableQ] Contour plots of the 1.3 cm and 7 mm continuum emis- 
sion are presented in Figs. [3^5] and [7] At 1.3 and 3.6 cm, the 
emission is resolved out by the VLA-A and we were not able to 
accurately calibrate the relative position of the continuum with 
respect to the H2O masers. On the other hand, the radio con- 
tinuum appears as a compact (or marginally resolved) source 
with the VLA-C and the peak positions of the 1.3 and 3.6 cm 
emissions coinc ide with one anothe r within the uncertainties. In 
agreement with Zapata et al. (2006), the peak of the 7 mm con- 
tinuum is found offset by l"-2" to the NW with r espect to the 
1.3 an d 3.6 cm emission peak. The 7 mm maps of Zapat aeTaT] 
(2006, Figs. 4 and 5), with a synthesized beam of about l'/7, 
show an unresolved spur of emission pointing to the SE toward 
the 1.3 cm continuum. We did not use their tapering of the data 
and thus attained at 7 mm a synthesized beam of about 077. 
With this resolution, the emission spur is resolved into a double 
lobe structure elongated in the SE-NW direction. This extended 
emission is detected at a low-significance level of o nly 3-4<x, but 
the fac t that is also observed in the tapered maps of lZapata et alj 
(2006) indicates that it is real. Figures [3^5] show that the 7 mm 
southeastern lobe is found at a position close to the peak of the 
1.3 cm (and 3.6 cm) continuum. We consider the small offset 
(between 071-072) of the 7 mm southeastern lobe with respect 
to the 1.3 cm continuum peak not to be real, but being an arti- 
fact owing to the low SNR of this 7 mm emission feature. In the 
following discussion, we assume a single source to be respon- 
sible for the emissions at 1.3 cm, 3.6 cm, and from the 7 mm 
southeastern lobe, and denote this source with the label "bl" (see 
Table [1]). Figure [2] shows the spectral energy distribution (SED) 
of the radio emission of source "bl". In TableQ] the 7 mm lobe, 
offsets to the NW of "bl" by about 075-077, is indicated with 
the label "b2". 



5.2. 22.2 GHz H 2 masers 

We searched for water maser emission through the whole field of 
view, l"xl ", and range of LSR ve locities, from 50 to 68 km s" 1 , 
explored bv lBeuther et al.1 (|2002c) using the VLA-B configura- 
tion. The centroid of their water maser distribution falls inside 
the area of our detections. 

We have detected 40 distinct water maser features, dis- 
tributed within an area of about 074 x 078 (Fig.[3j)). Most of the 
maser emission comes from 35 features distributed across a re- 
gion of size ~ 350 mas to the west of the 1.3 cm continuum peak. 
A cluster of 5 features spread over ~ 30 mas is observed offset 
by about 076 to the north of the 1.3 cm continuum source. The 
individual features properties are presented in Table[2] Maser in- 
tensities range from 0.10 to 65.49 Jy beam -1 . 14 features (35% 
of the total) persisted over at least 3 epochs, 9 of which lasted 4 
epochs. The spread in LSR velocities ranges from 67.9 km s _I , 
for the most redshifted feature (# 23), to 51.4 km s _1 , for the 
most blueshifted one (# 9). Water maser absolute positions and 
LSR velocities are plotted in Fig. 03. The inset of Fig. [3}? shows 
the details of a persistent structure of water maser emission, 
which presents a regular variation of LSR velocities with posi- 
tion, from 62 to 67 km s across a region of about 45 mas. This 
structure harbors several among the strongest maser features, 
including feature # 1 which experienced, between the first and 
second epoch, a powerful maser flare. Using a time baseline of 
9 months, we have measured relative transverse velocities of wa- 
ter maser features with a mean accuracy of about 30%. The am- 
plitude of relative proper motions ranges from 15.1+4.1 km s , 
for feature # 22, to 117. 1+7.0 km s _1 , for feature # 25. Figure^ 
presents relative positions and velocities of the features, deter- 
mined with respect to the bright, slowly-variable (< 10%) fea- 
ture # 4. Both the direction and the uncertainty of the relative 
proper motions are indicated. Figure|5j5 shows the derived abso- 
lute proper motions of the persistent water maser features. Since 
absolute proper motions are affected by large uncertainties, we 
show only the inferred mean direction of the motion. 

5.3. 1 .665 GHz OH masers 

iForster & Caswell (fl999h observed the HMSFR G16.59-0.05 
with the VLA (HPBW of 772 x 170) in the LCP band and re- 
ported 7 distinct centers of 1 .665 GHz OH maser emission. Their 
brightest detection had a flux density of 2.2 Jy at the Vlsr of 
58.9 km s _1 . Figur e |H shows the area of th e 1.665 GHz OH 
maser detections bv lForster & Caswell d!999t) . 

No 1 .665 GHz OH maser signal is detected on single VLBA 
baselines over 2 min scans, with an estimated rms sensitivity of 
0.9 Jy. The mean total-power spectrum of the 10 VLBA antennas 
(Fig. []]) shows a tentative detection in the RCP band of a 0.85 Jy 
line at a LSR velocity of about 59 km s _1 . By imaging maser 
data after phase-referencing to the calibrator J1825-1718, no 
signal was detected above a 5cr level of 0.06 Jy beam -1 , across 
an area of about 12" x 12" centered on the VLA 1.66 5 GHz OH 
maser position reported bv lForster & Caswell (QJ99). Since it is 
unlikely that the VLA OH maser position is wrong by more than 
6", we consider residual ionospheric phase errors as the most 
probable cause of image degradation and non-detection of the 
maser signal. 

5.4. 6. 7 GHz CH 3 OH masers 

We covered the range of LSR velocities (from 52 to 69 km s _1 ) 
over which 6.7 GHz maser emission was detected in the high 
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sensitivity, Parkes-64 m observations by ICaswell et al.l (Qj>95b), 
and imaged the whole fi eld of view (0723 X0 727) where CH 3 OH 
masers were detected by Wal sh et alj ([1998) using the ATCA in- 
terferometer. The centroid of the ATCA 6.7 GHz maser distribu- 
tion falls inside the area of our maser detections. 

We have detected 39 distinct 6.7 GHz CH3OH maser fea- 
tures, most (95%) of which distributed over an area of about 
074 x 074 centered on the 1.3 cm continuum peak. Hereafter 
we refer to this group of 6.7 GHz maser features as the "main 
cluster". Figure [3^ shows the spatial distribution and the LSR 
velocities of the 6.7 GHz masers. The individual feature proper- 
ties are presented in Table[3] Two features, labeled # 28 and 38, 
are found offset from the main cluster to the west and the 
north, respectively, and fall close (within about 30 mas) to water 
maser features with also a good correspondence in LSR veloc- 
ity (within 1.5 km s ; see Fig. @}. The spread in LSR veloc- 
ities ranges from 68.5 km s _1 , for the most redshifted feature 
(# 8), to 52.9 km s" 1 , for the most blueshifted one (# 38). The 
main cluster presents a bipolar distribution with a clear sepa- 
ration of l.o. s velocities: redshifted features toward the NW, and 
blueshifted ones toward the SE. The average l.o.s. velocities over 
the NW and SE regions are +3.7 and —1.5 km s _1 , respectively. 
Maser intensities range from 0.06 to 20.62 Jy beam -1 . 25 fea- 
tures (64% of the total) persisted over the 2 yr interval covered 
by our 3 observing epochs. For each maser feature, Table[3]re- 
ports the mean brightness variability defined as the ratio between 
the variation of the brightness (of the strongest spot) and the av- 
erage of the minimum and maximum brightness. Relative posi- 
tions are given with respect to the structurally-stable, compact, 
and bright feature # 2, whose intensity varied less than about 8% 
(see Fig. [6j». The mean brightness variability of methanol maser 
features was less than about 22%, spanning a range from 2%, for 
feature # 16, to about 45%, for feature # 21. These estimates of 
brightness variability should be taken as upper limits since they 
include effects from the slightly varying beam shape among dif- 
ferent epochs as well as amplitude calibration uncertainties. 

We have calculated the geometric center (hereafter "center 
of motion", identified with label # in Table [3]) of features with 
a stable spatial and spectral structure persisting over the 3 ob- 
serving epochs. We refer our measurement of transverse motions 
to this point (Fig. [7j>. Using a 2 yr time baseline, relative trans- 
verse velocities of individual maser features are determined with 
a mean accuracy of better than 30%. Proper motions, derived 
only for maser features with a stable (spatial and spectral) struc- 
ture, are calculated by performing a linear fit of feature position 
vs. time. An analysis of deviations from linear motion, to estab- 
lish the existence of possible accelerated motions, is postponed 
after an upcoming fourth epoch. The absolute distribution and 
relative proper motions of methanol 6.7 GHz masers are plotted 
in Fig. [7] The amplitude of relative proper motions ranges from 
3.6 ± 1 .4 km s -1 , for feature # 5, to 1 1 .2 ± 1 .6 km s -1 , for fea- 
ture # 6. The mean transverse velocity is equal to 6.2 km s , 
about 2-5 times greater than the spread in l.o.s. velocities. It is 
interesting to note that most of the transverse velocities for fea- 
tures placed in the NW of the main cluster are oriented towards 
N or NE, whereas features in the SE of the main cluster move 
towards S or SW. 

5.4.1 . Interpretation of 6.7 GHz maser proper motions 

In order to interpret correctly the proper motion of a maser fea- 
ture we need to establish if the measured motions are real or 
apparent. Change in physical and/or excitation conditions of 
maser-emitting gas can in principle lead to apparent motions 



(Christmas-tree effect). Concerning 22.2 GHz H2O masers, ac- 
curate proper motion measurements towards several HMSFRs 
strongly suggest that this maser emission traces physical mo- 
tions of gas bullets, po wered by stellar wi nds or jets emitted from 
massive YSO(s) (e.g.. lGoddietaT]|2006h . Since our study is one 
of the first attempts to use proper motions of 6.7 GHz methanol 
masers to infer gas kinematics in HMSFRs, we propose a few 
arguments supporting the kinematic interpretation of features' 
proper motions. 

Figure [6] shows the time evolution during 2 years (from 
February 2006 to March 2008) of the morphology and spectral 
structure of three persistent 6.7 GHz maser features (# 2, 8, and 
16), selected to be representative of the observed feature prop- 
erties. For each feature and each observing epoch, three sets of 
data are presented: the spectral profile (left panel of Fig. [6|; the 
image of the most extended feature emission obtained mapping 
the weakest channels in the spectral wings (hereafter the feature 
"image"; middle panel of Fig. |6)); the spatial and LSR velocity 
distribution of the brightest spots corresponding to the feature 
(hereafter the "internal velocity gradient"; right panel of Fig.®, 
To compare images of the same feature among different epochs, 
we cleaned the visibility data using the same restoring beam. 

The slightly resolved shape of the feature image remains ap- 
proximately constant across the epochs. That is essentially a con- 
sequence of the remarkable persistency of the feature spectrum 
and spot relative positions, with internal velocity gradients of 
typically ~ 0.1 km s _1 mas -1 (see Fig . |6j» . This might be inter- 
preted as the maser emission at different epochs coming from 
the same blob of gas, whose internal (physical and kinematical) 
properties do not change significantly over time. On the contrary, 
if proper motions were due to different regions of the cloud being 
excited at consecutive times, one would not expect the internal 
velocity and geometrical structure of the feature to be so well 
preserved. Forthcoming VLBI epochs will allow to establish the 
persistency of the feature spectrum and internal velocity gradient 
over a longer time baseline. 

Another indication that we are measuring real motions 
comes from the smooth variation of the observed 6.7 GHz maser 
proper motions with feature position. Looking at Fig.|7]one notes 
that proper motions of nearby features have similar orientations 
and amplitudes, which suggests that masers trace a smoothly 
varying velocity field. In contrast, the Christmas-tree effect is ex- 
pected to produce apparent motions uncorrelated with positions. 
Finally, considering the whole distribution of measured 6.7 GHz 
proper motions, it is quite evident that the velocities in the plane 
of the sky of redshifted and blueshifted features are antiparal- 
lel and perpendicular to the line connecting them (see Fig |7). 
In the next section, we interpret the whole distribution of l.o.s. 
and transverse velocities in terms of rotation around the center 
of motion of the 6.7 GHz masers. 

6. Discussion 

The simplest interpretation for the bipolar distribution of 
6.7 GHz maser l.o.s. velocities and the regular pattern of ob- 
served transverse velocities is that the masers are tracing rotation 
around the center of motion (see Fig. [7j». Since transverse veloc- 
ities are generally larger than l.o.s. velocities, most of the fea- 
tures move close to the plane of the sky. The bipolar distribution 
of l.o.s. velocities indicates that the rotation axis is inclined with 
respect to the line of sight, with redshifted (blueshifted) features 
to the NW (SE) of the main cluster rotating away from (towards) 
the observer. Taking the average inclination with the plane of the 
sky of all the measured proper motions, the rotation axis should 
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form an angle of about 30° with the l.o.s. . Assuming centrifugal 
equilibrium the required central mass is 35 M . This dynamical 
mass (Mdyn = RV 2 ot /G) has been obtained taking the mean ve- 
locity amplitude, V rot = 7.2 km s , and the average distance 
from the center of motion, R = 600 AU, of the 6.7 GHz maser 
features with measured proper motions. Although the derived 
value is affected by large errors owing to the uncertain source 
distance and the poor knowledge of the maser geometry, it is 
consistent with the 6.7 GHz masers being commonly found as- 
sociated with massive YSOs. The position angle -45°) of the 
main cluster of 6.7 GHz masers is approximately perpendicular 
to the NE-SW di rection of the l ess co llimated 12 CO (2-1) out- 
flow detected by iBeuther et all d2006l) . The NE-SW elongated 
spur visible in the NH3 (2,2) map of Fig.|4]might also mark dense 
gas close to the axis of this outflow. The elongated distribution 
of 6.7 GHz masers could then trace gas in a disk/toroid rotating 
around the massive YSO which drives the NE-SW 12 CO (2-1) 
outflow (sketch in Fig. [8]). 

Figure |7] shows the 6.7 GHz CH3OH masers lying on top of 
the 1.3 cm and 7 mm continuum source labeled "bl" in Table Q] 
The continuum spectrum (shown in Fig. [2]) is consistent with the 
source being an ultracompact (UC) H 11 region, with a turnover 
frequency of about 22 GHz, ionized by a Bl star. Note, however, 
that to fit the continuum spectrum a very small size of the UCH 11 
region is required, of the order of a few hundredths of arcsec 
(red dotted line in Fig.|2]i. That is in contrast with the continuum 
emission being not detected by our 1.3 and 3.6 cm VLA-A ob- 
servations (see Table [TJ- Since in our fit we have considered an 
homogeneous distribution of ionized gas, one possibility is that 
the continuum source is not a classical H 11 region, but that is 
characterized by strong density and (maybe) velocity gradients. 
We interpret the continuum source close to the 6.7 GHz masers 
in terms of free-free emission from the interaction of a jet, pow- 
ered by the YSO at the center of the 6.7 GHz maser distribution, 
with dense gas. Using the fluxes measured with the VLA-C at 
three wavelengths (3.6, 1.3, and 0.7 cm), the derived spectral 
index for the "bl" emission is a — 0.5 (see Fig. [2j, which is 
consistent with the interpretation in terms of a thermal jet. 

For shock induced ionizati on in a thermal jet, following the 
calculation by Anglada (1996) for optically thin emission, one 
has: Fd 2 = 10 35 (£2/47r) P, where F is the measured continuum 
flux in mJy, P is the jet momentum rate in M0 yr km s , Q 
is the jet solid angle in sr, and d is the source distance in kpc. 
Using the flux of 1 .4 mJy measured at 7 mm with the VLA-C 
for the source "bl", and the near kinematic distance of 4.4 kpc, 
one derives: P = 8 X 10~ 3 (Q/4^) _1 Mq yr~' km s _1 . The mo- 
mentum rate thus depends on the estimate of the jet collima- 
tion factor (Q/47r) _1 . Taking the center of motion of methanol 
masers as the YSO position, the (sky-projected) distance to the 
1.3 cm continuum source (about 071) and its size (> 1") im- 
ply an opening angle 9 > 2.0 arctan(0.5/0.1) = 2.7 rad, and 
a corresponding solid angle Q. > 5 sr. Thus, by assuming a 
wide-angle wind with Q. « An, the derived momentum rate is 
P = 8 x 10~ 3 Mq yr" 1 km s" 1 . From the value of the momen- 
tum measured by Beuth eret al.l ((2006, Table 5), using a dynam- 
ical time scale of ~ 10 4 yr (as derived from the outflow size of 
~ 10" and the average gas velocity of ~ 20 km s~' ), the esti- 
mate of the momentum rate for the NE-SW 12 CO (2-1) outflow 
is P = 6x 10~ 3 Mq yr" 1 km s" 1 . Such a value is consistent with 
that derived for the wind supposed to excite the source "bl", 
and supports our interpretation that the YSO associated with the 
6.7 GHz masers can be the one driving the NE-SW 12 CO (2-1) 
outflow (see Fig. [8]). 



The 7 mm VLA source "b2" is found offset by about 075- 
077 to the NW of the source "bl" (see Tableland Fig.0), in the 
region where emission of both 22 GHz water m asers and sev- 
eral h igh-density molecular tracers are observed dBeuther et al.l 
2006, Fig. 13). The continuum and molecular emission in this 
region might be powered by the same wind, emitted by the YSO 
associated with the 6.7 GHz masers, which is supposed to power 
the source "bl". However, a simple evaluation of the momen- 
tum rate necessary to drive water masers, leads us to exclude this 
possibility. Assuming that the observed water maser emission is 
excited in a jet from a star, knowing the average distance of wa- 
ter masers from the star and the average maser velocity, one can 
estimate the momentum rate of the jet (in the assumption that it 
is momentum-driven) with the expression: 

P = 1.5 X 10" 3 V 2 R 2 00 (Q./4tt) M Q yr 1 km s" 1 

where V10 is the average maser velocity in units of 10 km s" 1 , 
R100 is the average distance of water masers from the star in 
units of 100 AU, and Q. is the solid angle of the jet. In deriv- 
ing this expres sion, following mo dels of water maser excitation 
in fast shocks (Elit zur et alJ ll989). we have adopted an H2 pre- 
shock density of 10° crrT 3 . Taking the average distance of the 
detected water masers from the center of motion of the 6.7 GHz 
masers, corresponding to ss 600 mas (or 2600 AU), the av- 
erage (absolute) maser velocity of 60 km s" 1 , and Q. ^ An 
for a wide-angle wind, the derived momentum rate is P » 
36 Mq yr -1 km s" 1 . This value is exceedingly too high, about 
four orders of magnitudes higher than that derived above from 
the 7 mm emission of source "bl", and indicates that the YSO 
associated with the 6.7 GHz methanol masers is too distant from 
the water masers to drive their motion. 

We postulate the presence of another YSO driving the mo- 
tion of 22 GHz masers, powering the 7 mm continuum source 
"b2" and exciting the emission of the high-density molecu- 
lar tracers. It could be possibly located close to the cluster of 
strongest water masers (inset in Fig. [3), which presents the 
largest and most scattered proper motions (see Fig.|5]l. The or- 
dered distribution of l.o.s. velocities with position observed for 
this threadlike structure of water maser features, as well as then- 
fast variability (on a time scale between 3 to 6 months), lead us 
to interpret this emission in terms of a shock front driven by the 
nearby YSO. Both warm dust and ionized gas can contribute to 
the continuum source "b2", but having no information on the 
spectral index prevents us from distinguishing between the two 
contributions. The presence of strong water masers and the exci- 
tation of molecules as CH3CN, HCOOCH3, and HNCO, typical 
of HMCs, suggests that this YSO can be a massive one. 

The YSO exciting the 22 GHz masers is probably embed- 
ded in a denser cocoon of gas and dust than that harboring the 
6.7 GHz masers. The presence of water masers indeed witnesses 
high gas density with nn 2 ~ 10° cm 4 . The gas density is 
likely to increase along the SE-NW direct ion till reaching th e 
main mm continuum peak (source "a" of Zap ata et al. 2006), 
which marks the position of the most embed ded source in the re- 
gion. According to the model pro posed bv lBeuther & She pherd 
(2005) (see also lArce et al.ll2007D . outflows emitted by massive 
YSOs might increase their opening angle in the course of the 
YSO evolution, excavating larger and larger cavities in the YSO 
natal cocoon. Following this model, more embedded massive 
YSOs are likely to be also less evolved. In this view, one can 
speculate that the YSO at the center of the 6.7 GHz maser distri- 
bution (which we identify as the driving source of a poorly col- 
limated 12 CO (2-1) outflow observed at larger scale) is at a later 
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stage of evolution than that exciting the water masers, which 
in turn might be more evolved than the YSO associated with 
the main mm continuum emission. The most collimated (along 
the SE-NW direction) 12 CO (2-1) outflow of this regi on, re- 
vealed by the SMA observations of Beuther "et al.l d2006l) . could 
be driven by either the YSO close to the water masers or the 
one responsible for the main mm emission (sketch in Fig. [8]). 
Therefore, across a distance of * 0.05 pc, there are indications 
of sequential star-formation with younger and younger objects 
going from SE to NW. 

7. Summary and Conclusions 

Using the VLBA and the EVN interferometers, we observed the 
high-mass star-forming region G16.59-0.05 in the three most 
powerful maser transitions: 22.2 GHz H2O, 6.7 GHz CH3CH, 
and 1.665 GHz OH. The radio continuum emission toward 
G16.59-0.05 was also observed with the most extended VLA 
configuration to compare its brightness structure with the VLA- 
C archival data available. From our observations, we draw the 
following main conclusions: 

1. In the present work, we have collected evidence that 
methanol masers, alike water emission, are suitable tracers 
of the kinematics around massive young stellar objects. With 
three VLBI epochs over a time baseline of 2 yr, we have 
measured accurate (< 30%) proper motions of the 6.7 GHz 
masers. The l.o.s and sky-projected velocity distribution of 
methanol masers indicates that they are rotating with an aver- 
age velocity of about 7 km s _1 at distances of about 600 AU 
from a YSO. The inferred dynamical mass (35 M Q ) is con- 
sistent with 6.7 GHz masers tracing a high-mass YSO. The 
elongated NW-SE distribution of 6.7 GHz masers suggests 
they can originate in a flattened structure (disk/toroid) and 
that the YSO at the center of their distribution can be the 
one driving the motion of the NE-SW 12 CO (2-1) outflow 
observed at larger angular scale. 

2. Close to the 6.7 GHz masers, a compact (or slightly re- 
solved) continuum source is observed with the VLA-C at 
3.6 and 1.3 cm, and at 7 mm. We interpret this radio contin- 
uum in terms of free-free emission of gas ionized by a wide- 
angle wind emitted by the YSO associated with the 6.7 GHz 
masers and interacting with the surrounding dense gas. The 
momentum rate of this wind is consistent with that measured 
for the NE-SW 12 CO (2-1) outflow. 

3. Water masers are found offset by more than Of! 5 from the 
center of the 6.7 GHz maser distribution, in the same re- 
gion where a weak 7 mm VLA source and emission in high- 
density molecular lines are observed. The derived absolute 
velocities of the water masers, although affected by large er- 
rors, seem to indicate fast motions with an average amplitude 
of about 60 km s . We postulate the presence in this region 
of a distinct YSO, responsible for driving the motion of wa- 
ter masers and exciting the continuum and molecular line 
emissions. 

Previous SMA observations in several molecular lines have 
revealed multiple molecular outflows from the G16.59-0.05 
star-forming region, suggesting the presence of several massive 
YSOs across a region of « 0.1 pc. This work shows that maser 
VLBI data are a powerful tool to identify the massive YSOs driv- 
ing the large scale outflows, to study the gas kinematics close 
to the forming stars, and to detect disks/toroids of size of hun- 
dreds AU at typical distances of several kiloparsecs. In the next 



future, when VLBI maser observations will be complemented 
with data of new generation (sub)millimeter interferometers (i.e. 
ALMA) at comparable angular resolution, it will be possible to 
better constrain the physical and kinematical properties of both 
outflow(s) and core(s) hosting maser activity. 
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Fig. 1. Total-power spectra of the H2O, CH3OH and OH masers toward G16.59-0.05. Upper panel: system-temperature (T sys ) 
weighted average of the 22.2 GHz total-power spectra of the 9 VLBA antennas observing on 2007 January 4 (BR, HN, KP, LA, 
MK, NL, OV, PT, SC). Middle panel: Effelsberg total-power spectrum of the 6.7 GHz methanol maser emission on 2007 March 16. 
Lower panel: T SJJ -weighted average of the 1.665 GHz total-power spectra of the 10 VLBA antennas observing on 2007 April 13. 
Continuous (black) and dashed (red) lines are used to denote the right (RCP) and left (LCP) circular polarizations, respectively. The 
vertical dotted line across the spectra indicates the systemic velocity (V sys ) inferred from CH3CN measurements. 
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Table 1. G16.59-0.05: Radio continuum associated with the CH3OH and H2O maser emission. 
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(1) Observations from archival data, project code AZ149 (P.I. Zapata). 

(2) Data from the CORNISH survey: "The CORNISH Survey of the Galactic Plane", C. R. Purcell, M. G. Hoare, P. Diamond, 2008, Astronomical 
Society of the Pacific Conference Series, Vol. 387, Pg. 389. 
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Fig. 2. Spectral energy distribution of the VLA component "bl" toward the HMSFR G16.59-0.05. Dots and errorbars report the 
values and the associated errors (l<x) of the fluxes at 0.7, 1.3 and 3.6 cm measured with the VLA-C (see Table [TJ. The derived 
spectral index, indicated by the black dotted line, is a = 0.5. The red dotted line shows the best fit of the measured fluxes with a 
model of an homogeneous H n region ionized by a Lyman continuum of 3 x 10 45 s _1 with a radius of about 0V04. 
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Fig. 3. Absolute positions and LSR velocities of maser species observed in G16.59-0.05: a) CH3OH (dots); b) H2O (triangles). 
Different colors are used to indicate the maser LSR velocities, according to the color scale on the right-hand side of the plot, with 
green representing the systemic velocity of the HMC. The VLA 1.3 cm and 7 mm continuum emissions are plotted with dotted and 
dashed contours, respectively. The 1.3 cm contour levels range from 40 to 90% of peak emission (0.6 mJy beam -1 ) at multiples of 
10%, and from 70 to 90% for the 7 mm peak (0.7 mJy beam -1 ). Inset: enlargement of the N-S threadlike structure described in 
Sect. 15.21 Maser features are indicated with the label numbers given in Table [2] 
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Fig. 4 . Collection of the sub-arcsec observations toward G16.59-0.05. Left panel: map in the NH3 (2,2) line from Codella et al. 
(119971) . Contour levels range from 30 to 90% of the peak emission (38 mJy beam ') at multiples of 10%. The restor ing beam is 
shown in the lower left corner of the panel. The star marks the peak position of t he 7 mm continuum emission imaged by Zapat a et al.l 
( 2006). The empty square indicates the area where OH masers were detected by Fo rster & Caswelll (Il999l) . Right panel: enlargement 
of the region over which we have detected maser emission. Red triangles and blue dots represent H2O and CH3OH maser positions 
from our VLBI measurements, respectively. The VLA 1.3 cm and 7 mm continuum emissions are given with dotted and dashed 
contours, respectively, using the same contour levels shown in Fig. [3] 



14 



A. 



Sanna et al.: Studying Massive Star Formation through maser VLBI 





Fig. 5. 22.2 GHz H2O maser kinematics toward G16.59-0.05. a) positions (triangles) and transverse velocities of the H2O maser 
features relative to the feature #4 (indicated with the vertex-connected symbol). Colored cones are used to show both the direction 
and the uncertainty (cone aperture) of the proper motion of maser features. The proper motion amplitude scale is given by the 
black arrow on the bottom left corner of the panel. Different colors are used to indicate the maser LSR velocities, according to the 
color scale on the right-hand side of the panel, with green denoting the systemic velocity of the HMC. The VLA 1.3 cm and 7 mm 
continuum emissions are given with dotted and dashed contours, respectively, using the same contour levels shown in Fig. [3] b) 
absolute positions and transverse velocities of the water maser features. The plotted field of view is the same as in the upper panel 
and symbols and contours have the same meaning as in the upper panel. Absolute transverse velocities of water masers are affected 
by large uncertainties and only the mean direction (and amplitude) of motion is shown. 




Fig. 6. Upper, middle, and lower set of plots present the time evolution of the spatial and spectral structure of three selected 6.7 GHz CH3OH features. Left to right plots present 
the spectral profile, the image, and the internal velocity gradient of the feature, respectively, at our 3 observing epochs (corresponding to the upper, middle, and lower panels 
of each set of plots). Spectral profile: dots report the intensities of feature's spots, emitting at different Vlsr; vertical dashed lines delimit the velocity ranges of the spectral 
wings, over which maser signal has been summed to produce the feature image shown in the middle plot. Image: map of the spectral-wings emission of the feature, obtained by 
summing channel maps over the two velocity ranges indicated in the spectral profile. The origin (0, 0) of each map is centered on the peak position of the brightest spot of the 
feature. Contour levels range from 30 to 90%, at multiples of 10%, of the peak intensity reported on the top right corner of each panel. At each epoch, the same restoring beam 
has been used and is shown in the lower left corner of the lower panel. The maps in the lower set of plots show the emission of three closeby, spatially blended, features, indicated 
with the label numbers reported in Table [3] To recover the most extended structure of this spatially-blended emission, in this case we have produced tapered maps ("UVTAPER 
= 10 MA"). Internal velocity gradient: spatial and LSR velocity distribution of the brightest spots of the feature. Different colors are used to indicate the maser LSR velocities, 
according to the color scale on the right-hand side of the plot. Symbol size scales logarithmically with spot intensity. In each panel, spot positions are measured from a point u, 
(0, 0) which coincides with the origin of the feature image for the corresponding epoch. 



16 



A. Sanna et al.: Studying Massive Star Formation through maser VLBI 




L^j — i ^ , , 1 , ^ , d 1 68 9 

18 21 m 09 s 14 09 s 12 09 s 10 

RIGHT ASCENSION (J2000) 

Fig. 7. 6.7 GHz CH3OH maser kinematics toward G16.59-0.05. Absolute positions (dots) and transverse velocities of the CH3OH 
maser features relative to the center of motion (as defined in Sect. \5A\ of the methanol maser distribution (indicated by the cross). 
Colored cones are used to show both the direction and the uncertainty (cone aperture) of the proper motion of maser features. The 
proper motion amplitude scale is given by the black arrow on the bottom right corner of the panel. Different colors are used to 
indicate the maser LSR velocities, according to the color scale on the right-hand side of the plot, with green denoting the systemic 
velocity of the HMC. The VLA 1.3 cm and 7 mm continuum emissions are given with dotted and dashed contours, respectively, 
using the same contour levels shown in Fig. [3] Numbers close to proper-motion cones of a few features are the feature labels reported 
in Tableland mark the features whose properties are presented in Fig. [6] 
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Fig. 8. Schematic picture of the main emission components detected toward the HMSFR G16.59-0.05, within a field of view of 
about 10". The drawing is not to scale. The sources labeled "a", "b " (splitted in t his pa per in two components, "bl" and "b2"), and 
"c" mark the three VL A radio continuu m components detected by Zapa ta etaD {2006). The big grey circles identify the two mid- 
IR so urces detected bv lDe Buizer et all (120051) . The north western one is a l so associated with the hot ammonia core (|Cpdella et al.l 
1 19971) and molecular, thermal and maser, line emissions (IWalsh et al.lll998HForster & Caswellfl 999: Beut her et all2 006). The main 
directions of the two CO bipolar outflows (NE-SW and NW-SE) resolved bv lBeuther et al.l (120061) are shown. The distribution and 
main direction of motion of the methanol (blue and red dots according with approaching and receding l.o.s. velocities) and water 
(black dots) maser components reported in this paper are presented: 1) the CH3OH masing-gas traces a disk/toroid rotating about a 
massive YSO, powering the NE-SW outflow and a thermal jet associated with the continuum source "bl"; 2) the H2O masing-gas 
together with the radio component "b2" and the thermal lines traces a distinct massive YSO, possibly in an earlier evolutionary 
phase than the one exciting the methanol masers. 
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Table 2. Parameters of VLBA 22.2 GHz water maser features. For each identified feature, the label (given in Col. 1) increases with decreas- 
ing brightness. Cols. 2 and 3 report the LSR velocity and brightness of the brightest spot (for the epoch specified by the number in brackets). 
Cols. 4 and 5 give the position offset relative to the feature # 4, toward the East and the North directions, respectively. Cols. 6 and 7 report the 
projected components of the feature proper motion relative the the feature # 4, along the East and North directions, respectively. Col. 8 shows the 
relative error of the proper motion amplitude. The absolute position of the reference feature # 4 is reported at the bottom of the Table. 
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1 oa cc , n no 

120. jj ± u.uy 








Zj 


AA T/1 
00. /4 


r\ a o /o\ 
U.4y (2) 


1 1 c o/i _i_ n no 
— 1 15.84 ± U.U8 


1 ct io _i_ n 1 o 
lj /.32 ± U.12 


OT A _i_ c n 

— 2/.4 ± j.U 


111 O _i_ T A 

1 13.8 ± i.o 


AO/ 

Oto 


ZO 


A/1 1 f\ 
04. 1U 


u.40 (I; 


1 nA a n _i_ n 1 n 
— 1U0.4 / + U. 1U 


1 it c/i _i_ n 1 o 
13 /.j4 ± U.12 








in 

Z 1 


AA 10 
00.32 


r\ /I c /o\ 


1 n z;i , n no 
— 11 /.0j + U.U8 


1 AO A O _i_ n 1 o 
102.42 ± U. 12 








1Q 

Zo 


AC A 

Oj.48 


r\ /] r» /o\ 
U.4U (2J 


1 m nA _i_ n no 

lu/.uo ± u.uy 


1 no /i n _i_ n in 
iuy.4u ± U.IU 








Zy 


AO HA 

02. /4 


u. jy i-j j 


1 1 o nn _i_ n no 

— i i2.uu + u.uy 


iioic_i_ni/i 
1 18. 3j ± U.14 








30 


60.30 


0.35(1) 


6.92 ± 0.09 


12.23 ±0.12 


4.9 ±3.6 


20 .4 ± 4.6 


21% 




An Al 

0U.03 


H 'XA ( A\ 

UJ4 (4J 


lU/.lO ± U.l 1 


111 ^A xh 16 

1 1 1.j4 ± U.lo 








32 


62.20 


0.33(1) 


-109.48 + 0.19 


128.05 ± 0.20 








33 


60.84 


0.25 (4) 


5.81 ±0.09 


10.10 ±0.14 








34 


60.84 


0.21 (3) 


-8.41 ± 0.09 


13.03 ±0.20 








35 


62.32 


0.21 (4) 


-111.81 ±0.08 


120.87 ±0.12 








36 


63.46 


0.18(1) 


-107.52 ±0.10 


131.67 ±0.13 








37 


60.09 


0.16(1) 


106.91 ±0.10 


119.47 ±0.15 








38 


58.61 


0.10(1) 


228.34 ±0.10 


-31.43 ±0.18 


12.8 ±4.3 


10.4 ±7.0 


34% 


39 


63.04 


0.10(1) 


264.41 ±0.11 


684.83 ±0.15 








40 


60.00 


0.10(3) 


108.53 ± 0.09 


113.03 ±0.15 














Reference Feature: Absolute position 








Feature 




R.A.(J2000) 
(h m s) 




Dec.(J2000) 

C ") 







4 



18:21:09.0921 1±0.00009 



-14:31.48.6790±0.0023 
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Table 3. Parameters of EVN 6.7 GHz methanol maser features. For each identified feature, the label (given in Col. 1) increases with decreasing 
brightness. Cols. 2, 3 and 4 report the LSR velocity and brightness of the brightest spot, and its percent brightness variability. Cols. 5 and 6 give the 
position offset relative to the feature # 2, toward the East and the North directions, respectively. Cols. 7 and 8 report the projected components of 
the feature proper motion relative to the center of motion (as defined in Sect. 15.41 and identified with label # 0), along the East and North directions, 
respectively. Col. 9 shows the relative error of the proper motion amplitude. The absolute position of the reference feature # 2 is reported at the 
bottom of the Table. 



Feature 


^LSR 


Fpeak 


Var. fl 


Ax 


Ay 




Vy 


A|V|/|V| 




(kms -1 ) 


(Jy beam -1 ) 




(mas) 


(mas) 


(km s" 1 ) 


(km s" 1 ) 




1 


62.11 


20.62 


20% 


-106.27 + 0.09 


205.81 ±0.13 








2 


59.47 


12.90 


8% 








-3.2 ± 1.0 


-4.0 + 1.4 


25% 


3 


64.13 


7.88 


20% 


-72.93 + 0.09 


175.07 ± 0.13 








4 


58.42 


7.60 


3% 


40.35 ± 0.09 


-18.67 ± 0.13 


-1.7 ± 1.0 


-6.1 ± 1.4 


22% 


5 


63.78 


7.15 


35% 


-158.52 + 0.09 


135.34 + 0.13 


-0.6+ 1.0 


-3.5 + 1.4 


40% 


6 


66.06 


4.15 


9% 


-206.96 ± 0.09 


113.62 ± 0.13 


1.8 + 1.1 


11.0+ 1.5 


14% 


7 


67.81 


2.98 


18% 


-200.44 ± 0.08 


120.97 + 0.12 


0.7 + 1.0 


4.1 + 1.4 


35% 


8 


68.52 


2.42 


18% 


-198.32 + 0.10 


120.56 ± 0.14 


5.4+ 1.2 


3.3 + 1.7 


22% 


9 


63.95 


2.03 


17% 


-58.07 ± 0.10 


177.76 ± 0.17 








10 


59.91 


1.77 


45% 


37.02 + 0.10 


3.00 + 0.15 


-4.2 ± 1.2 


-6.1 + 1.8 


22% 


11 


64.74 


1.60 


32% 


-212.91 + 0.10 


125.86 ± 0.16 








12 


57.02 


1.20 


21% 


44.19 + 0.09 


-24.65 ± 0.16 


-1.3 ± 1.2 


-4.6 ± 2.0 


41% 


13 


58.42 


1.16 


8% 


20.18 + 0.18 


-18.57 ± 0.46 








14 


64.30 


1.13 


30% 


-85.95 + 0.10 


168.74 + 0.17 








15 


57.45 


0.99 


10% 


-100.07 + 0.10 


-52.87 ± 0.16 








16 


56.84 


0.87 


2% 


-82.83 + 0.11 


-55.85 + 0.18 


-6.3 ± 1.6 


-2.4 ± 2.6 


27% 


17 


62.55 


0.72 


44% 


-210.41 + 0.13 


101.17 ± 0.27 








18 


58.51 


0.56 




-15.95 + 0.22 


0.44 ± 0.33 








19 


58.60 


0.49 




69.45 + 0.19 


-77.32 + 0.31 








20 


61.84 


0.41 


25% 


-102.92 + 0.15 


189.05 ± 0.30 








21 


61.76 


0.40 


45% 


-81.04 + 0.19 


112.36 ± 0.46 








22 


58.86 


0.28 


34% 


-6.48 ± 0.18 


-30.55 ± 0.39 








23 


64.04 


0.18 




-68.46 ± 0.42 


103.02 ± 0.58 








24 


59.12 


0.17 




84.43 ± 0.39 


-78.86 ± 0.66 








25 


60.09 


0.17 


10% 


-11.76 + 0.22 


-25.28 ± 0.50 








26 


58.68 


0.16 




-17.03 ± 0.63 


-94.58 ± 0.91 








27 


59.12 


0.16 


17% 


26.70 ± 0.24 


-18.15 ± 0.63 








28 


58.07 


0.14 


19% 


-396.02 + 0.25 


-39.64 ± 0.64 


... 






29 


63.07 


0.13 


44% 


-213.83 + 0.28 


195.95 ± 0.63 


... 






30 


57.63 


0.12 




-62.37 ± 0.30 


-57.29 ± 0.58 








31 


61.76 


0.12 




-68.25 ± 0.42 


93.82 ± 0.77 








32 


60.70 


0.11 




-287.88 + 0.56 


-43.17 + 0.80 








33 


58.77 


0.09 




-7.14 + 0.97 


68.86+ 1.10 








34 


61.32 


0.09 


15% 


-68.81+0.33 


122.23 ± 0.80 








35 


63.25 


0.08 




-157.22 + 0.65 


33.89+ 1.01 








36 


63.51 


0.08 




-209.00 + 0.54 


104.94+ 1.03 








37 


63.69 


0.08 




-149.29 + 0.55 


90.60 + 0.57 








38 


52.89 


0.07 




-407.11+0.65 


669.44+ 1.07 








39 


56.84 


0.06 




-91.01+0.95 


-124.54+ 1.31 

















-84.71+0.06 


67.90 + 0.10 











Reference Feature: Absolute position 



Feature 


R.A.(J2000) 


Dec.(J2000) 




(h m s) 


C ") 


2 


18:21:09.13711+0.00023 


-14:31:48.6853±0.0036 



(a) Variability (%) = (F M ax - F M in of the brightest spot) / ((F M ax + F M in)/2.0) 



